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We present an experimental realization of an autonomous Maxwell’s Demon, which extracts microscopic 
information from a System and reduces its entropy by applying feedback. It is based on two capacitively coupled 
single electron devices, both integrated on the same electronic circuit. This setup allows a detailed analysis 
of the thermodynamics of both the Demon and the System as well as their mutual information exchange. The 
operation of the Demon is directly observed as a temperature drop in the System. We also observe a simultaneous 
temperature rise in the Demon arising from the thermodynamic cost of generating the mutual information. 


Thermodynamic processes are governed by fundamental 
laws, of which the first, conservation of energy, is paramount 
in all fields of physics and cannot be violated at any level of 
description known to date. The Second Law in turn states that 
entropy, the measure of disorder, of a closed system cannot 
decrease. This has most important consequences, such that 
heat flows from hot to cold, irreversible processes must dissi¬ 
pate work, and devices of perpetual motion are impossible. To 
challenge this law, James Clerk Maxwell presented a thought 
experiment in 1867 of a “finite being” capable of accurately 
measuring the velocity of molecules IT) . It would act between 
two separated reservoirs, permitting only fast molecules to en¬ 
ter one reservoir, while allowing only the slow ones to the 
other. Under such a process heat is transferred from cold to 
hot, apparently in violation of the Second Law. This idea, 
coined as ’’Maxwell’s demon” by Lord Kelvin, has over a cen¬ 
tury spurred further research on the relation between informa¬ 
tion and energy establishing quantitative relations ll2l 4T2ll . On¬ 
going progress in nanotechnology has also provided concrete 
means to test such relations experimentally GMl, thus re¬ 
igniting acute interest in actually constructing a Demon. 

Recently, several theoretical proposals on configurations in¬ 
cluding both the System as well as the Demon have been pre¬ 
sented ElElMl. Such a configuration is known as an au¬ 
tonomous Maxwell’s demon, for the fact that the measure¬ 
ment and feedback operation takes place internally. Here, we 
experimentally realize an all-in-one Maxwell’s demon, whose 
operation principle is cartooned in Fig.[T](a). The System is a 
single electron transistor (SET) lIZTl . formed by a small nor¬ 
mal metallic island connected to two normal metallic leads 
by tunnel junctions. The two junctions permit electron trans¬ 
port by tunneling, and are assumed to be identical (both with 
the same resistance Rf). The Demon measures the number 
of electrons on the System island, and applies feedback as 
depicted in Fig. [T] (a). When an electron tunnels to the is¬ 
land, the Demon traps it with a positive charge (panels 1 and 
2). Conversely, when an electron leaves the island, the De¬ 
mon applies a negative charge to repel further electrons that 
would enter the island (panels 3 and 4). The System elec¬ 
trodes contain a reservoir of conduction electrons whose ther¬ 


mal excitations provide sufficiently high energy carriers to 
overcome the the trapping or repulsion induced by the De¬ 
mon, contributing heat Q = —AE where AE is the energy 
cost of the tunneling event. In doing so, the System entropy 
decreases as ASg = Q/Tg, where Tg is the System reservoir 
temperature, i.e. the Demon extracts information of tunnel¬ 
ing electrons to apply feedback that causes the entropy of the 
System to decrease. While the configuration resembles the¬ 
oretical proposals on quantum dots 1^ 1^ 1^ . and shares 
features with the Coulomb drag effect ITOllTn . it constitutes a 
genuine autonomous Maxwell’s Demon where only informa¬ 
tion, not heat, is directly exchanged between the System and 
the Demon. 

Our experimental, autonomous realization of the cycle in 
Fig.[T](a) relies on coupling the System island capacitively to 
a single electron box (SEE), a small normal metallic island 
connected by a tunnel junction with resistance i?d to a single 
normal metallic lead. Here, the SEE undertakes the role of the 
Demon. The resulting Hamiltonian is 

H{n,N) = Eg{n-ngf+Ed{N-Ngf+2J{n-ng){N-Ng), 

( 1 ) 

where the dynamic variables n and N are the net number 
of electrons that have entered the System and the Demon is¬ 
lands, respectively. Eg and E^ are their charging energies, 
while J > 0 describes their mutual Coulomb interaction and 
is essential for the device operation. The state (n, N) evolves 
when an electron tunnels through a junction, n changes to 
n (T) 1 when an electron tunnels from (to) the System is¬ 
land. Correspondingly, N changes to N (T) 1 when it trans¬ 
fers from (to) the Demon island. Constant external control 
parameters Ug and Ng govern the System current and the cou¬ 
pling of the Demon to the System, respectively. The System 
is voltage V biased, such that the electron (with elementary 
charge — e) tunneling in the direction of (against) the voltage 
bias experiences an energy cost AE = AH (T) eV 12, where 
AH = iT(n ± 1, N) — H{n, N) for changing n is given by 
Eq. 0. Similarly, for the electron tunneling in the Demon 
AE = H{n, N±l)- H{n, N). 

The interaction between the System and the Demon is max- 
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FIG. 1: Operation principle, (a) The Demon monitors the System (a 
single-electron transistor) for electrons that tunnel into (panel 1.) or 
out of (panel 3.) the island. It then immediately performs a feedback 
by applying a positive charge to trap (panel 2.), or negative charge 
to repel (panel 4.) the electrons. Coulomb blockade ensures that 
only either one or zero electrons reside in the System island. The 
electrons are always tunneling against the potential induced by the 
Demon, and therefore the System cools down, (b) Energetics of the 
System under voltage bias V in the experimental and autonomous 
realization of the cycle in (a) with another single electron structure 
operating as the Demon. The conduction electrons of the System 
follow Fermi distribution, providing electrons that can overcome the 
energy cost J — eV 12, where J is the coupling energy between the 
System and the Demon, however in doing so the System cools down 
by an equal amount. The energy J is dissipated by the Demon as it 
reacts, changing the projected energy cost experienced by the elec¬ 
tron tunneling in the System from — J — eV/2 to J — eV 12. Note 
that here the System island is drawn without Fermi distribution for 
simplicity. Also, the described operation could in princple be per¬ 
formed non-autonomously by externally measuring the System state, 
and changing the energy of the System as feedback, see e.g. 1251 . 


imized by setting Ug = Ng = 0.5, producing the Hamiltonian 
H{n,N) = J{2n — 1){2N — l)/2 and energetics depicted 
in Fig. [^(b). We furthermore require eVjksT Es,E^, 
such that only the lowest energy states of Eq. 0 are avail¬ 
able, such that both n and N are practically limited to two 
possible values, 0 and 1. States (n = 0, = 1) and 

(n = 1,7V = 0) are charge neutral, both with energy —J/2. 
Here, we refer to either of the states as ’ground’ or g. The 
state {n = 0,N = 0) has an overall positive charge and 
(n = 1, = 1) an overall negative charge. We refer to 

them as ’charged’ or c, both with energy J/2. Any single 
tunneling event will take p to c or c to g, with respective 
AHg^c = J = —AHc^g. We assume that the System is at 
uniform temperature Tg while the temperature of the Demon 
is Td, such that the occupation probability distribution Pn,N 
obeys Po,i = ^* 1,0 = Pg/‘^ and Pq.o = = -Pc /2 with 

Pg ~ ^ c—(r g^cP rc—^ 5 ) and Pq — rg->.c/(rg_>.c+rc_>.g). 
Here, with notation J± = J ± eE/2, the term Fc-^g = 
rs(—J_|_) + rs(—J_) + Td{—J) is the overall transition rate 
from c to g, while Fg-^c = Ts{J-) + PsiJ+) + ^d{J) is the 


corresponding overall transition rate from g to c as a sum of 
rates in the System in the direction of bias, against the bias, 
and the transition rate in the Demon, respectively. The transi¬ 
tion rates are 


Ts/di^E) 


1 AE 

e'^Rs/d - 1 ■ 


( 2 ) 


The charge current in the System is / = (e/2)(rs(J_) — 
rs{J+))Pg + (e/2)(rg(-J+) - Tsi-J-))Pc and the total 
heat generation rate there is 


Qs=- (J-Ts (J-) + J+Ps iJ+))Pg 

+ {J.Ts (-J-) + J+Fg i-J+))Pc, 


(3) 


reflecting the fact that if the Demon successfully maintains a 
high Pg by feedback as in Fig. [T](b), Qs is negative. Similarly, 
the rate of heat generation in the Demon is 

Qd = -JrdiJ)Pg + JTdi-J)Pc, (4) 


which in turn is positive as the Demon applies feedback on 
states c as in Fig. [T] (b). Consider Tg = Td = T. It 
can be shown that when feeTtanh (J/2fcBr) < (<^/4)(l + 
Rd/Rs)~^, Eq- 0 gives negative Qs, i-S- cooling, within a 
range of 0 < \V\ < |14iax| < 2 J/e (see Supplementary mate¬ 
rial for derivation). The entropy of the System then decreases 
as Ss = Qs/T < 0 seemingly against the Second law, how¬ 
ever we still get Sd = Qd/T > —Qs/T = Ss resulting from 
Joule’s law, Qs + Qd = IV. 

Although energetically our device follows Joule’s law, it 
is the information flow between the System and the Demon 
that permits the decrease of System entropy. The mutual 
information between the System and the Demon is = 
ln(P„,Ar) - ln(P„) — In(PAr), where Pn and TV are the oc¬ 
cupation probabilities of n and N, respectively. As Pn=o = 
Pn=i — Pn=o = Pn=i = 0.5, mutual information changes 
in a tunneling event from p to c as A/i„,g->c = ln(T’c/T’g), 
and for c —> g as A/m c->g = — A/m,g->.c Tun¬ 

neling events in the Demon change mutual information at the 
rate 

im,d = In ^3 + In Td i-J) Pm (5) 

Majority of the tunneling events in the Demon are c —>■ g 
transitions, and since Pg > Pc, Im,d is positive. The rate of 
mutual information change by the System tunneling events is 
Im,s = —Im,d- As discussed in Ref. Il33l . the System heat 
generation satisfies Qs > —fcBTg/m.d implying that the max¬ 
imum amount of cooling is bound by the amount of mutual 
information generated by the Demon. Correspondingly, gen¬ 
erating mutual information has a thermodynamic cost for the 
Demon as Qd > ^Bld/m.d- This can also be understood in 
terms of the configurational entropy 5'conf = — In (P(n, N)) 
as follows ll^ : tunneling events in the Demon bring the cir¬ 
cuit from unlikely state c to the more probable state g, decreas¬ 
ing 5'cont ■ At least an equivalent of heat must be dissipated to 
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FIG. 2: Experimental realization, (a) A scanning electron micro¬ 
graph of the structure. False colour identifies the System island (light 
blue), its left lead (dark blue), and right lead (dark green), as well as 
the Demon island (orange) and its leads (red). The System temper¬ 
ature deviations from their base value, ATl, ATr, and ATd, are 
measured at the indicated locations (see Supplementary Material for 
details of measurement setup), (b) I aiV = 120 ^V. When Ng is an 
integer, I is modulated by Ug as in a standard SET. When Ng ~ 0.5, 
I is smaller due to Demon interaction, (c) ATd at V = 120 ^V. 
When Ug, Ng ~ 0.5, ATd elevates due to the information flow be¬ 
tween the System and the Demon. Measured data in (b, c) are shown 
on the left and numerically obtained predictions on the right. 


satisfy the second law. On the other hand most of the tunnel¬ 
ing events in the System bring the setup to a more improbable 
state c, increasing conhgurational entropy. The Second law 
then allows cooling by at most the amount of configurational 
entropy decreased, i.e. —ASs < ASconf- We note that in 
the limit <C Rs, P{n,N) follows the thermal equilibrium 
distribution of the Demon. Then \n{Pg/Pc) = J/ksT^ such 
that Ijn.d = Qd/feeTd by Eqs. Q and ([^. This implies that 
measurement of heat generated in the Demon is also a direct 
measurement of information extracted by the Demon. 

Figure (a) shows a scanning electron micrograph of the 
experimental realization of the Maxwell’s demon. It was fab¬ 
ricated by standard electron beam lithography combined with 
shadow evaporation 0^ of copper (normal metal) and alu¬ 
minum (superconductor) metal Aims. Our device has the fol- 
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FIG. 3: Operation as a Maxwell’s demon and as a one-sided refrig¬ 
erator. Quantities shown are I (black), ATl (blue), ATr (green), 
ATd (red), with parameter values F = 20 pV, ro,s = 77 mK, and 
ro,d = 55 mK. (a) Measurement at Ng = 0.5 (Maxwell’s demon). 
Both Tl and Tr decrease, indicating overall cooling of the System. 
This is justified by the mutual information transfer between the Sys¬ 
tem and the Demon, which in turn generates heat in the Demon, ob¬ 
served as elevated Td. (b) Measurement at Ng — 0 (SET refrig¬ 
eration (21 [351). Either ATl or ATr can be negative, however, 
not simultaneously: overall heat is generated in the System. Mea¬ 
sured data (symbols) are shown on the left and numerically obtained 
predictions (lines) on the right, (c) Energetics at different operation 
points, indicated as numbers in panels (a) and (b). At the operation 
point 1, the Demon is interacting with the System as in Fig. 1 (b). At 
operation points 2-4, the Demon is inactive. 


lowing parameters: Es/k^ — 1.7 K, — 810 mK, 

J/^B — 350 mK, Rs ~ 580 kQ, and i?d — 43 kfl (two 
parallel junctions each with ~ 85 kQ tunneling resistance). 
The fully normal System and Demon junctions are realized 
with laterally proximized aluminum dot technique OtI . We 
determine the heat generated in the left (L) and right (R) lead 
of the System as well as the lead of the Demon by measur¬ 
ing the respective temperatures Tl, Tr, and Td, as indicated 
in Fig. 1^ (a). This is achieved by reading the voltage of 
current-biased normal metal - insulator - superconductor junc¬ 
tions, see e.g. Ref. Ii38]l . Finally, the leads of the System and 
the Demon are interrupted with direct contacts to supercon¬ 
ducting leads, which permit charge transport by Andreev pro¬ 
cesses li39l but block heat transport at low temperatures. The 
structure is measured in a ^He / ^He dilution refrigerator at the 
bath temperature of 40 mK. Details on the device fabrication 
and measurement configuration are given in the Supplemen¬ 
tary Material. 

The continuous heat generation is mediated primarily by 
lattice phonons that couple with the conduction electron 
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heat bath at temperature TL/R/d, contributing Qmph = 
- Tl), m = L, R or d, where E is a material 
specific constant, Vm is the volume of the circuit element, and 
Tq „i is the base temperature Ho). For the left and right elec¬ 
trodes of the System, Vl/r ~ 2.8 /rm x 70 nm x 20 nm. Its 
island is approximately twice as large in volume. The Demon 
has the total volume Vd « 4 x 3.2 /im x 150 nm x 20 nm. 
We use E « 4 X 10® Wm“^K“® for Cu. The rate of electron 
tunneling (10® Hz) in our device is faster than the phonon re¬ 
laxation rate (10^ Hz), however it is small compared to the 
inelastic electron-electron relaxation rate, which is typically 
of the order of 10® Hz Bdl . allowing the electrodes to equi- 
libriate to an effective electron temperature that deviates 
from To m. Furthermore, the temperature change caused by 
an individual tunneling electron is sufficiently small so that 
QlTm is a good approximation for the entropy change. The 
temperature Tm equilibrates such that the net heat generation 
is zero, i.e. Qm = —Qph,m- The base temperature To.m is 
measured at Ug = Ng = 0, where the state is Coulomb block¬ 
aded to n = N = 0 corresponding to the energy minimum 
in Eq. 0 and no heat is generated in the circuit. Figure |^(b) 
shows that charge current / in the System modulates with Ug 
as in a standard SET. However, when Ng = 0.5, the maxi¬ 
mum measured current is reduced due to the feedback by the 
Demon. Eigure[^(c) demonstrates how at Ug = Ng = 0.5, 
the heat generated in the Demon is maximized for extracting 
information of the transported electrons. 

The main result of this paper is presented in Eigure (a), 
showing our observation at V = 20 p,V ~ 23jie of how 
the System cools down and its entropy decreases. Simulta¬ 
neously, we observe how the Demon, which collects the in¬ 
formation and immediately applies a feedback to the System, 
generates heat as a necessary thermodynamic cost for extract¬ 
ing information from the System. On the other hand, Eigure 
[^(b) shows unchanged at Ng = 0 since the Demon is ef¬ 
fectively uncoupled from the System as its state is locked to 
N = 0. With that Coulomb blockade refrigeration ll34l lT5l 
occurs when Ug deviates from 0.5 by causing either the left 
or right lead to cool down, but overall heat is generated and 
entropy is produced in the System. 

Eigure (a) shows a measurement of current (inset) and 
temperatures as a function of V at Ug = Ng — 0.5. Increas¬ 
ing voltage bias boosts electrons to pass through the System, 
however, at the cost of lower entropy decrease per electron. 
Eurthermore, the risk of electrons to pass through the System 
without feedback control from the Demon increases, in partic¬ 
ular via multi-electron tunneling (see Supplementary material 
for details). Eigure (b) compares the heat and mutual in¬ 
formation produced by the Demon, demonstrating that they 
differ by less than 15% for low V. The data shown in Eig. 
(c) shows improvement of entropy decrease to up to 20 p,V, 
beyond which errors in the feedback process overcome the 
benefit of enhanced rate of electron injection. At this voltage, 
the cooling power on the System is estimated to be —Qs ~ 6 
aW, while the heat dissipation in the Demon is Qd « 19 aW. 
Based on the heat generation, the mutual information produc- 



FIG. 4: Bias dependence. Here, rig = 0.5, Ng = 0.5 and To.d = 55 
mK. The data points (symbols) are obtained by averaging over 210 
repetitions, (a) ATl (blue squares), ATr (green circles), and ATd 
(red diamonds) with their respective prediction with To,a = 77 mK 
(dashed lines) and Tq.s = 62 mK (solid lines). Inset: I in the same 
measurement. Applying voltage increases the number of electrons 
passing through the System and in turn the information flow between 
the System and the Demon. This is observed as increased Td. (b) 
Numerical comparison between Qd/ksT^ and Im,d, demonstrating 
that the two quantities match, (c) Enlarged view of the measured 
ATl (blue squares) and ATr (green circles). Increasing voltage 
bias further enhances the entropy decrease in the System to up to 
about ±20 fjN. The model assumes a perfectly symmetric System 
and therefore predicts equal ATr and ATr with the fit To^g = 62 
mK (solid line). 

tion rate by the Demon is then /i„,d « 25 x 10® Hz. The 
current is / « 600 fA, i.e. ^ 4 x 10® electrons cross the 
System per second. Should successful feedback be performed 
for every electron, the heat extracted by the Demon would be 
I X 23/e Ki 36 aW. Experimentally we extract « 52%, of this 
value i.e. this fraction of the electrons transported through the 
System are successfully feedback-controlled by the Demon. 
Eor efficiency at maximal cooling power, —Qs/IV, we the 
get Ri 0.56. 

In conclusion, we have realized and demonstrated experi¬ 
mentally a physically transparent autonomous Maxwell’s De¬ 
mon on a chip, based on coupled single-electron circuits un¬ 
dergoing tunneling events in a self-controlled manner. The 
Demon acts on the System to decrease its entropy, observed 
as a temperature drop. The configuration allows one to mea¬ 
sure the effect of the Demon on the System, as well as to mea¬ 
sure the thermodynamics of the Demon itself. The device pre¬ 
sented here demonstrates how information is transferred from 
the System to the Demon, leading to heat generation in the 
Demon in amount that corresponds to the rate of information 
transfer. This setup constitutes a step towards autonomous 
information-powered nanodevices. 
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Supplementary Material 


Fabrication and measurement setup: 



FIG. SI: Device fabrication and measurement configuration. False colour identifies the metal layers evaporated during the fabrication. Light 
blue depicts Al, and light orange is the Cu evaporated prior to oxygen exposure. Any contacts between these two layers are fully transparent. 
The final red is the Cu evaporated after the oxygen exposure. Any contacts to Al are tunnel junctions, whereas contacts to the first Cu layer are 
still highly transparent due to its slow oxidation rate. Fgi and Vg 2 determine the control parameters rig and Ng. By applying 7th ,l ~ 7th ,r ~ 
7th,d ~ 1 pA, Kh,L, Vth,R, and Lth,D depend approximately linearly on the temperatures Tl, Tr, and Td. A 7heat,L ~ 7heat,R ~ 10 pA 
current is applied over the heater leads to bring Tr and Tr to the temperature of the System island. 

The System and Demon were formed simultaneously in four steps. First, a 15 nm Al (superconductor, S) layer was evaporated 
to form the bottommost shadow, see Fig. |S1| This was directly followed by an evaporation of 20 nm Cu (normal metal, N) layer 
to form the middle shadow. Any metal-to-metal contacts between Al and Cu are transparent NS contacts. After the second 
evaporation, the sample was exposed to oxygen, forming a thin AlxOy layer on any Al not already covered by Cu. Finally, a 
second 20 nm layer of Cu was evaporated as the topmost shadow. Most of the contacts formed between Al and Cu at the third 
evaporation are NIS tunnel junctions with the AlxOy as the insulator (I). However, small Al dots with a direct contact to Cu 
are rendered normal due to inverse proximity effect llJTlI . This technique is used to form the NIN junctions for the System and 
the Demon. Copper oxidizes significantly slower than Al, thus any Cu-Cu contacts are practically transparent. The capacitive 
coupling between the islands of the SET and the Demon was achieved with parallel plate capacitance to an underlying Au strip 
evaporated prior to the SET - detector structure, covered with « 30 nm of AlxOy achieved by atomic layer deposition (ALD). 
As the NIS junctions are operated as thermometers, they simultaneously induce some cooling ll3^ on the electrodes probed. To 
initialize the System in such a way that Tr and Tr are at the same temperature as the island, the System leads are equipped with 
additional current-biased tunnel junctions that induce heating to compensate for this NIS cooling effect. 


Tbnneling and heat generation rates 

This section gives the equations for the tunneling rates through the junctions in this device. Tunneling events may take place 
through the left or right System junction both with a resistance 7?s ~ 580 kH, or through the left or right Demon junction both 
with resistance 7? ^ 85 kH. However, as the Demon is not voltage biased, the two junctions are parallel junctions to the same 
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potential, and can thus can be modeled as a single one with i?d ~ 43 kfl. A tunneling event may take place between the left 
System lead (L) or right System lead (R) and the System island (I) at rate 

r„j(Ai;) = ^ J */m(e)(l -Me- AE)), (6) 

where AE is the energy cost for the tunneling electron, m = L, R, and fk{e) = (1 + exp(e/fcB7fe))~^ is the Fermi-Dirac 
distribution of the conduction electrons in the metal fc = L, R, I. Similarly, the tunneling rate between the Demon lead (d) and 
its island (D) is 

^ J de/d(e)(l - Mie - AE)). (7) 

We assume that Td = M, in which case Eq. (0 simplifies to rd(A£') = rd,D(Ai7) = (e^i?d)~^— 1). 
Moreover, if the System is at uniform temperature implying Tl = Ti = Tr = Tg, the System tunneling rate simplifies to 
T,{AE) = Ti^.iiAE) = Tnj{AE) = {e^R,)-^AE /_ i). 

Next follow the equations for the rates of heat generation. Heat generation rate in the metal to = L, R by electrons tunneling 
between to and I is g™ i(Ai7) = {e^Rs)~^ f de(—e)fm(e)(l — Me — AE)), and in I it is i(Ai7) = (e^i?s)“^ / de{e — 
AE)fMe)il - Me - AE)). With Tr = Ti = Tr = Tg, these simplify to q^jiAE) = qij{AE) = q^J^E) = qkj(AE) = 
— ( 2 e^Rs)~^AE^/(e^-^ABTs _ = q^(AE)/2, where MAE) is the total rate of heat generation on the System, summing 

over the rates of the lead and the island between which the tunneling takes place. It also satisfies 

MAE) = -AEMAE). (8) 

Finally, the heat generation rate by electron tunneling in the Demon (summing over d and D, with Td = Tr) is qd{AE) = 
-(e2i?d)-iA£;7(e^-®/'=B'rd - 1) = -ATrd(AT). 

We note that the equations above only depend on the energy cost for the tunneling process. It depends on the initial state 
(n, N) and source and destination of the electron, as 

^ 2 [Tg(n - ng) + J(iV - A,)] ^ , 

^^R^i(i^R) ^ 2 _ n,) + J{N - NM M ^ 

^^d^D(D^d) ^ 2 [Td(A - M) + J{n - n,)]. 


Full expressions for I and Q 


Although the states {n,N) = (0, 0), (0,1), (1,0), and (1,1) are dominant, we also consider higher energy states in the 
simulations. A master equation for the configuration at steady state is written as 


dPn.N 

dt 


- [M,{AE]M^) + rL.i(ATi>7 + Tr^atR^I) + rR,i(AT^7f) + rd(ATd7D) + rd(A£;D77]p„,^ 

+ [Mi{AE]MIn)^^K,i{AE^M,n)\ Pn-1,N + [M,l{AEi-^lM +MA^Ei-^fM] Pn+l.N 

+ rd(AT^7° + rd(ATD77)P„,jv+i = 0, (10) 


which can be solved from Eq. ( [T0| ) by noting n En,N = 1- In the steady state, the current through the System 
is the same in each cross-section and can then be written as / = [rL,i(AT');'7) — rL,i(AT77)] Pn,N- The 

rate of heat generation in the metal to = L, R is Qm = Pn.N, and in I it is 

Qi = Y.n,N + 7l.i) + 4,1 + 4,i(^^4^)] Pn,N, and finally in the Demon it is Qd = 

Sn N + 7d(^T'n 7'^)] Pn,N■ The simplified forms written in Eqs. (2, 3) in the main text are obtained in the 

two-by-two-state limit n, A = 0,1 at = Ag = 0.5, with uniform temperatures Tg and Td, Qs = Ql + Qi + Qr, using Eq. 

The final temperatures are determined by solving the heat balance equation for each circuit element, indexed by to = L, R, I, 
or d, as Qm{Tm) + Qm,ph{Tm) + Qm,heater{Tm) = 0. The term Qm,heater(Tm) is Otherwise zero but for the metal to = L, 
R, Qm.heater = ~ ^ is the heat leak through the heater junctions by Wiedemann-Franz law, Lq « 2.44 x 10“® 

WHK“^ is the Lorenz number, T^ is the electrode temperature, To^m is the base temperature of that electrode, and Ph ~ 2 MH 
is the resistance of the heater junctions. 
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Temperature threshold for cooling 


In this section we derive the threshold temperature for the Demon to be able to cool down the System. We consider Ug = 
Ng = 0.5 and uniform temperatures Tg for the System and for the Demon. In the main text Eq. (2), the rate of heat generation 
in the System is 

Qs = -J-Ts (J-) Pg - J+Eg (J+) Pg + J+Eg (-J+) Pe + J-Eg i-J-)Pc, (11) 

where J± = J ± eE/2. With notation = e^PgEg(a:) = x/ (exp(^sx) — 1), = e^PsEd(a;), and /3s = l/k^Ts, the second 

derivative of Eq. •HD is 


dPQ. 

dV^ 


v=o 


2RsJ{2gj + 2g_j + Dj + D_j) = Agjg_jJj3s - D_jgj [2(1 - (isd-jT “ JPld-j] + 


Dj9-j [2(1 - PsdjY - + 4^sffjg-j (-2 + /3gpj + /3sP_j) 


gjD-j - g-jDj 

‘^g-j + D-j + 2gj + Dj 


If we further have Tg = Td = T, we can write Dj = {Rs/Rd)gj and Eq. ([T^ reduces to 


d^Qs 

dV^ 


Psgjg-j 


y=o 2Rs{2 + Rs/ Rd){gjg-j) 


4+^(4-J/3gC0th(4/3gJ)) 


( 12 ) 


( 13 ) 


t 2 I 

The critical threshold to achieve cooling is ^ 0’ i-®- ,3s<^coth (/3g J/2) > 4 (1 + Pd/Ps); in which case for a finite 

V we can achieve negative Qg. In the limit Pd ^ Rs, the numerically obtained threshold is k^T < 0.2611 x J. 


Multi-electron tunneling 

Tunneling beyond single-electron processes is considered in the same way as in Ref 041 . In particular, we consider the 
multi-electron tunneling for the data in Eig. 4 of the manuscript, where Ug = Ng = 0.5. When n = 0 and iV = 1 (or n = 1 
and N — 0), any single electron tunneling event costs energy. Eor example, a tunneling event in the System from source to 
the island would bring n from 0 to 1 with an energy cost J — eV 12. However, if another electron would afterwards tunnel in 
the Demon, bringing N from 1 to 0, that electron would gain J in energy, resulting in a total energy gain eV 12. Eurthermore, 
in our experiment the Pd — 43 kH is close to the resistance quantum, Pk ~ 26 kH, giving rise to multi-electron tunneling in 
scenarios, where single electron tunneling costs energy whereas with multiple tunneling events the total energy cost would be 
negative. Consider an electron tunneling event with an energy cost AP. Right after the event, there are six possible follow-up 
transitions (indexed fc), 1: L —> /, 2: / —> P, 3: P —4: / —> P, 5: d —D, or 6 : D —d, each with its respective energy 
cost APfc (see Eq. ([^). In this context, the equations for calculating rates of tunneling and heat induced are modified by energy 
broadening 7 (e) Efc(APfc — AP — e), where Efe(APfe — AP — e) the tunneling rate of the transition k (see Eqs. to 

lzl)- The rate of heat generation is modified as 

gbr(AP) = i J deq{AE + e)^^^^^^, (14) 

where the subscripts and superscripts are the same for q^r and q on the left and right hand side of the equation (for example, 
to solve I br’ would insert i to the right-hand side). Similarly, the modified tunneling rates are given by Ebr(AP) = 
7 r“t J fier(AP + e) 7 (e)/(e^ + 7 (e)^). The rates evaluated by the this expression are used in Eq. ( [TOl i to evaluate the probability 
distribution Pn,N- Energy conservation demands that the total rate of heat generation ( 7 tot,br(AP) induced by the tunneling 
processes Ebr(AP) satisfies ( 7 tot,br(AP) = —AP Ebr(AP). In addition to the heat generation of Eq. ( [T4| , each junction k 
dissipates energy e as heat at a rate proportional to |Efe (AP^ — AP — e)/(e^ -f 7 (e)^). Thus the additional rate of heat generated 
by the virtual processes for each junction (including the junction where the actual tunneling event occurs) is virt(^£') = 
7 r“^ f deeE(AP -f e)Efc(AP/j — AP — e)/(e^ + 'y(e)^), approximately to split evenly between the electrodes shared by the 
junction. Now gtot,br(AP) = -f (jb° + 'Zfc.virt(AP) = — AP Ebr(AP), where ’from’ and ’to’ refer to the source and 

destination of the tunneling electron, satisfies energy conservation. Einally, we note that in the limit of Pp —00 with p = s, 
d, the effect of multi-electron tunneling should vanish: indeed with 7 (e) 0, (jbr(AP) —>■ q{AE) and Ebr(AP) E(AP), 

whereas 4 ^virt(AP) 0. 

In addition to energy broadening, we consider direct co-tunneling processes, where two electrons tunnel simultaneously. 
The first tunneling event takes place between pair (L,I), (R,I), or (d,D), which we indicate with A, and the second tunneling 
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FIG. S2: Multi-electron tunneling. Parameters are the same as in Fig. 4 of the main manuscript. Solid lines show predictions with multi¬ 
electron tunneling, dashed ones include only single-electron tunneling, (a) Cooling effect as a function of V. (b) Current through the System 
as a function of V. Multi-electron tunneling causes additional current to leak through the System, reducing the cooling power, (c) Temperature 
in the Demon as a function of V. 


event takes place over another pair, indicated with B. The energy cost for the events are AEa and AEb, respectively, in the 
case only that event took place. If both tunneling events take place, the total energy cost is AE. When AE < AEa + 


AEb, cotunneling events become relevant. The rate of co-tunneling is Tct = f AE+ ) 


heat dissipation rate in A is (7 a, ct = ^ / 


4A(|AB+e)rB(iAB-e) 


rA(|Ag-r£)rB(iAg-£) 
;a- 5 Ab-£)(A£;b- 5 A.e- 

, where the superscript for q is the same for left and 


{AEA-^AE-e){AE-B-^AE+e) 

right hand side of the equation (see the similar discussion after Eq. ([T^). The heat dissipation rate in junction B is qB,ct = 




rAi^AE+(.)qB(^AE-e) 


Figure 


S2 


shows the effect of multi-electron tunneling in contrast to what would be 


'(A£;a-5A£;-£)(A£)b-5A.E+£) ■ 

expected if only single-electron tunneling would take place. Additional current leaks through the System without the feedback 
from the Demon, resulting in smaller cooling power on the System. 


Gate dependence of temperatures 


Figure]^ shows how our configuration behaves as a function of Ug and Ng. Figure 3 in the main text has been obtained from 
these type of measurements. 
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FIG. S3: Dependence of control parameters Ug and Ng at V = 20pV for current I (first top and bottom panels, starting from left), ATl 
( second panels), ATr (third panels), and ATd (last panels). Measured values are shown on top panels and simulated data are shown at the 
bottom. It is possible to identify the regime where the Demon begins to interact with the System (identified by the crossing point of the dashed 
lines), observed as suppressed maximum I, vanishing of heat generation in the System, and heat generation arising in the Demon. 
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Discussion on demon memory 

Many theoretical works on autonomous Maxwell’s demons iCHni consider systems in contact with information reservoirs, 
represented by sequence of bits in the models. In these models the information reservoirs exchange no energy with the system. 

The relation between our setup and these models can be understood by considering the two-level demon state as a memory 
(one bit). In the other works the sequence of bits has been initialized to a preferred initial state. That is to say, the system 
has been subjected to a pre-determined stream of bits. In our setup the sequence of bits is prepared by the demon itself while 
it operates, and not given by initialization of the stream. Because the demon has to erase the single-bit memory for the next 
feedback operation, there is a thermodynamic cost for this erasure, which is seen as dissipation of heat in the setup. 

As the authors in Ref im point out, it is possible to study the entropic interactions between the system and the information 
reservoir using entropic currents, which is presented by the flow of mutual information in our model. In the works ll^ fTOll . the 
authors show in detail that in fact the entropy production and the dynamics can be understood either by considering the flow 
of mutual information or the model of information resource with entropy difference in the incoming and outcoming bit streams 
(entropy difference in the information resource). 

The change in the demon state (flipping of a bit) changes the energy barriers between the system states. At this stage, the 
system does not yet experience transitions and therefore does not cool down, thus there is no energy transfer between the system 
and the demon. The carriers of heat in the setup are the electrons, and no electrons can tunnel between the system and the demon, 
such that direct heat exchange is in fact impossible in our setup. 
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